Charge Order in NaV 2 5 studied by EPR 
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We present angular dependent EPR measurements in NaV20s at X-band frequencies in the tem- 
perature range 4.2 K < T < 670 K. A detailed analysis in terms of the antisymmetric Dzyaloshinski- 
Moriya and the anisotropic exchange interactions yields the following scheme of charge order: 
On decreasing temperature a quarter-filled ladder with strong charge disproportions, existing for 
T > 100 K, is followed by zig-zag charge-order fluctuations which become long-range and static 
below T SP = 34 K. 



The observation of an exponential decrease in mag- 
netic susceptibility in a '-NaV 2 5 below 34 K g] trig- 
gered many experimental and theoretical investigations. 
Following the determination of the crystal structure by 
Carpy and Galy j^j as space group V2\mmn with linear 
chains of V 4+ ions (spin S = 1/2) separated by non 
magnetic V 5+ chains, NaV^Os, analogous to CuGe03, 
was classified as inorganic spin-Peierls system. Later on 
a re-investigation of the crystal structure showed that 
NaV2C>5 has to be considered a quarter-filled ladder 
system |^-|| with only one vanadium site in the high- 
temperature phase described by space group Pmmn. The 
direction of the ladders is given by the crystallographic 
b axis. Each rung consists of two VO5 pyramids along 
the a axis, which share one corner at their base. Based 
on this structure a charge-order transition followed by 
different kinds of spin order were proposed [ppp]. A 
recent investigation suggested a low-temperature phase 
consisting of modulated ladders with zig-zag charge order 
alternating with non modulated ladders with vanadium 
ions of intermediate valence V 4 5+ ||. 

Electron-paramagnetic resonance (EPR) is ideally 
suited for the investigation of the magnetic properties of 
vanadium systems, since the vanadium ions themselves 
(as V 4+ with electron configuration 3d 1 in the present 
case) can be used as a microscopic probe of the spin sys- 
tem. In this paper we present angular dependent EPR 
measurements at X-band frequencies (9.48 GHz) within 
the temperature range 4.2 K < T < 670 K. The experi- 
mental details concerning the crystal growth of NaV^Os 
and EPR measurements have been described in previous 
papers || . For all temperatures the EPR spectrum con- 
sists of a single strongly exchange-narrowed Lorentzian 
line. The EPR intensity does not show any angular de- 
pendence within experimental error. Its temperature be- 



havior follows that of a one dimensional antiferromag- 
netic Heisenberg chain at high temperatures T > 100 K 
fll(| , whereas it decreases exponentially below the transi- 
tion temperature Tgp. Strong deviations from the model 
predictions were observed between 100 K and Tgp ||. 
The g-values are found as g a ~ <?b ~ 1-98 and g c f=a 1.94, 
which are typical for V 4+ in octahedral crystal symmetry 
|jll|l , and increase only slightly below Tsp- Here we con- 
fine ourselves to the discussion of the angular and temper- 
ature dependence of the resonance linewidth, which car- 
ries the information about the interactions of the vana- 
dium spins with their local environment. 

Figure [l] shows the temperature dependence of the 
linewidth AH with the external magnetic field applied 
along the three main axis of the crystal. Starting with 
an isotropic value AH = 8 Oe at the transition temper- 
ature Tsp = 34 K, the linewidth increases monotonously 
with increasing temperature and develops a remarkable 
anisotropy with respect to the c axis of the crystal. In the 
respective temperature regime the curvature of AH{T) 
starts with a positive sign and changes to negative val- 
ues above 100 K. Towards low temperatures T < Tgp 
the line broadens anisotropically again probably due to 
the unresolved hyperfine structure or due to a spin-glass 
transition pl| . 

It has been reported by Yamada et al. |l3| that the line 
broadening in the high-temperature regime T > 100 K 
must be due to the antisymmetric Dzyaloshinski-Moriya 
(DM) interaction, which is the only possible mechanism 
to explain the observed order of magnitude of some 
100 Oe for the linewidth and its anisotropy. Estimations 
for anisotropic exchange and dipole-dipole interactions 
yield values which are 100 and 1000 times smaller, re- 
spectively. The angular dependence of the linewidth is 
shown in the upper inset of figure |l| at four different tem- 
peratures, where the b axis was taken as rotation axis 



1 



perpendicular to the static magnetic field. 
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FIG. 1. EPR linewidth AH as function of temperature T 
for the magnetic field applied along the three crystallographic 
axes. Upper inset: Angular dependence of the EPR linewidth 
AH for the magnetic field applied perpendicular to the b axis 
at temperatures T > 100 K. Lower inset: Temperature depen- 
dence of the strength parameter Adm of the DM interaction 
(see equation |l| and text). 



Following Yamada the data are well described by 
AH = A DM • (1 + cos 2 1?) + AH a 



(1) 



where •& is the polar angle with respect to the c axis. Adm 
is proportional to the strength of the DM interaction 
and AHq the residual linewidth due to further relaxation 
mechanisms. Above 100 K the parameter Adm strongly 
increases, as it is shown in the lower inset of figure |l| 
and the experimental ratios AH C / AH a and AH C / AH^ 
nearly approximate the value of 2, underlining the dom- 
inant influence of the DM interaction. 

Below T = 100 K the results strongly deviate from the 
high-temperature behavior. The positive curvature of 
the linewidth indicates a change of the dominant inter- 
action. In the temperature range 34 K < T < 60 K, the 
maximum linewidth is not found for the magnetic field 
applied in c direction any more, but it appears with re- 
spect to the a axis. Figure nicely shows this cross-over 
regime. The angular dependence within the b-c plane 
shown in the upper inset of figure || is of special interest, 
because here an additional modulation appears, which is 
well described by 



The amplitude A of the modulation is about 1 Oe, which 
is of the order of magnitude estimated for the anisotropic 
exchange interaction, and increases slightly with decreas- 
ing temperature, whereas the prefactor B, which can 
be considered to be due to the DM interaction, using 
cos(2$) = 2 cos 2 ?? — 1, as discussed below, strongly de- 
creases. 

The angular dependence of the EPR linewidth contains 
important information about the local electronic distri- 
bution on the vanadium ladders. We start our discussion 
with respect to the high-temperature regime T > 100 K, 
where the linewidth is determined by the DM interaction 



Hi 



dij ■ [Si x Sj 



(3) 



The vanadium spins Si and Sj are coupled by superex- 
change via an oxygen ion. The direction of the DM vector 
dij is determined by 



dy = d ■ [riio x n Q j], 



(4) 



where the space vectors nio and noj connect the spins i 
and j with the oxygen-bridge ion respectively. The maxi- 
mum linewidth, which assigns the direction of the DM 
vector, is found for the magnetic field applied parallel to 
the crystallographic c axis. 
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FIG. 2. Cross-over regime in the temperature dependence 
of the EPR linewidth AH. Upper inset: Angular dependence 
of the linewidth AH for the magnetic field applied perpen- 
dicular to the o axis at temperatures T < 100 K. Lower inset: 
Fit parameters A and B derived from equation ^| as function 
of temperature T. 



AH = Acos{Ad) + Bcos(2d) + AH 



(2) 
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FIG. 3. Projection of the ladder structure of VO5 pyramids 
into the a-b plane for different possible charge distributions, 
the spheres denote the vanadium ions, the oxygen ions have 
been omitted for simplicity. The oxygen ions of the relevant 
V-O-V bridges are marked by the DM vectors dij. (a) Sym- 
metric charge distribution (grey clouds) of a quarter-filled spin 
ladder V 4 5+ ; (b) chains of V 4+ (dark) and V 5+ (light); (c) lo- 
cal charge disproportion at high temperatures; (d) charge or- 
der fluctuations (dynamic, short-range) above Tsp and charge 
order (static, long-range) below Tsp. Case (c) and (d) are in 
full agreement with the EPR results. 



Therefore the relevant oxygen-bridge vectors riio and 
noj must build an angle smaller than 180° within the 
a-b plane. Moreover the existence of a non vanishing 
DM interaction requires a unit cell without inversion cen- 
ter, because otherwise the DM vectors cancel each other, 
as it is the case for a quarter-filled spin ladder, where 
each vanadium spin is equally distributed between two 
V 5+ ions on the rungs of the ladder (cf. Fig. ||a). The 
structure determined by Carpy and Galy |2j fulfills the 
condition of asymmetry, because of the separation of the 
vanadium ladders in V 4+ and V 5+ chains. Here the DM 
interaction takes place only via one oxygen bridge be- 
tween two neighboring V 4+ sites along the b axis (cf. Fig. 



||b). However, using the structural data ||| for the vana- 
dium and oxygen ions under consideration, we find that, 
with respect to the vanadium positions, the oxygen ions 
within the chains are farther displaced along the c axis 
Ac w 0.58 A than along the a axis Aa « 0.28 A. This 
yields a dominant contribution of the DM vector point- 
ing into the a direction, which is not observed experimen- 
tally. We conclude that the EPR data strictly contradict 
both proposed high-temperature structures: the linear 
spin chain j| as well as the symmetric quarter-filled lad- 
der HH . 

Recently, detailed optical investigations by Damascelli 
et al. |l4|] revealed a local charge disproportion on each 
rung of the ladders, which also destroys the inversion 
symmetry yielding a finite DM interaction. Hence the 
charge distribution has to be taken into account more 
carefully (cf. Fig. p^) . Following the analysis of Dama- 
scelli, the center of the electronic distribution x c on the 
rungs of the ladder can be determined. Considering each 
rung as an independent linear molecule of two V 5+ ions, 
the eigenstates of the additional electron in between are 
given by = uif)\ + vip T , where the electron stays near 
the left-hand side and ip2 = uip r — vipi, where the elec- 
tron favors the right-hand side. The atomic orbitals ip T 
and ipi °f the right-hand and left-hand side are occupied 
according to the following weights: 
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(5) 



The parameter A describes the asymmetry of the onsite 
energy of the electron located in i/> r with respect to 
The charge-transfer energy Ect = VA 2 + At 2 is deter- 
mined by the hopping integral t and resembles the energy 
splitting between the two eigenstates. We choose the ori- 
gin of the coordinate system x = in the left-hand ion, 
giving the position of the right-hand ion at x — R. If the 
electron favors for example the left-hand side, the center 
of the charge distribution x c can be obtained from the 
eigenstate ipi as 



v 2 R + {R/2)2uvS Rh 
u 2 + v 2 + 2uvSrl 



2 Ect 



) (6) 



where Srl is the overlap integral between ip T and ip\ . The 
last expression was obtained assuming «Sr,l ^ E which is 
justified with respect to the distance R « 3.6 A of both 
ions, and using equation |^. Taking the experimental val- 
ues A w 0.8 eV and t ss 0.3 eV from Ej, we calculate 
x c = 0.1R rj 0.36 A, yielding a shift of the electronic 
charge distribution by 0.36 A with respect to the cen- 
ter of the left-side V 5+ ion. If now neighboring rungs of 
the ladder locally obey the same charge distribution, we 
obtain an overall shift of Aa « 0.64 A of the vanadium 
spins with respect to the connecting oxygen ion. This 
yields an V-O-V bridge angle of about 140° within the 
a-b plane and therefore a dominant DM contribution in 
c direction. The EPR data provide strong experimental 
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evidence for the appearance of charge disproportions in 
the quarter-filled spin ladder at high temperatures. 

The peculiarities observed in the cross-over regime 
(34 K < T < 60 K) can be described in the framework of 
the competition of antisymmetric DM exchange and sym- 
metric anisotropic exchange, which has been discussed by 
Yamada et al. pq| for the one dimensional antiferromag- 
net KCUF3 and more extensively for ferromagnetic Cu 
layers by Soos et al. Q . In these papers the evolution of 
the linewidth with temperature is calculated for both in- 
teractions. In both cases the linewidth saturates at high 
temperatures: The antisymmmctric DM exchange yields 
a strong temperature dependence of the linewidth with a 
negative curvature within the whole temperature regime. 
The anisotropic exchange gives rise only to a weak tem- 
perature dependence of the linewidth with a positive cur- 
vature at low temperatures changing to a negative one 
at high temperatures. Combining both contributions, we 
qualitatively obtain the temperature dependence of the 
linewidth observed in NaV^Os experimentally. Compar- 
ing these data to those of the one dimensional Heisenberg 
antiferromagnet KCuF 3 in more detail, we recognize an 
important difference. Whereas in KCUF3, in agreement 
with theoretical predictions, the DM interaction vanishes 
exactly at T = 0, for NaV2 05 the extrapolation of the 
high-temperature data suggests the DM interaction to 
disappear already at 100 K, as documented in the lower 
inset of figure [p. This indicates that already for temper- 
atures T < 100 K the charge disproportions change to 
a symmetric zig-zag order of V + and V 5+ ions in the 
ladder, where the DM vectors cancel each other (cf. Fig. 
|d). This result is in good agreement with reports of 
charge-order fluctuations for 80 K > T > Tgp from Ra- 
man spectroscopy Jjj]] and the onset of a static long-range 
zig-zag order below Tsp for ultrasonic jL8j and dielectric 
results H). 

Finally the 7r/2-periodic modulation of the angular de- 
pendence in the b-c plane according to equation @ can be 
understood in terms of the anisotropic exchange, only. 
The DM interaction produces 7r-periodic modulations 
in any case. The anisotropic-exchange interaction yields 
7r/2-periodic modulations only if its secular contributions 
are enhanced with respect to the non-secular contribu- 
tions by spin diffusion, which is characteristic for low 
dimensional systems 16 . As we will report in a more 



detailed subsequent publication, the fact that the tt/2- 
periodic modulation is only observed in the b-c plane, 
whereas it does not occur in both a-b and a-c plane, 
means that all elements of the anisotropic exchange ten- 
sor vanish except Jb c 7^ 0. Under these conditions, the 
prefactors of equation || resemble anisotropic A ~ j£ c 
and antisymmetric B ~ ^4dm exchange. From the lower 
inset of figure @ we obtain that for 60 K > T > 34 K 
the Dzyaloshinsky-Moriya interaction more and more 
decreases with decreasing temperature whereas the 
anisotropic exchange slightly increases. The anisotropic 
exchange parameter Jb c is due to the coupling between 



electron spins of adjacent layers in the c direction. This 
is a strong hint that the ordering below Tgp involves 
vanadium ions of neighboring layers. 

In conclusion the EPR data confirm the structure 
of a quarter-filled spin ladder with strong local charge 
disproportions at high temperatures T > 100 K, which 
result in a non- vanishing DM interaction (Fig. |^c). The 
direction of the DM vector is determined by the V-O- 
V bridge angle ss 140° along the chains within the a-b 
plane. The pronounced weakening of the DM interac- 
tion below 100 K but far above Tsp = 34 K indicates the 
onset of zig-zag charge-order fluctuations, which become 
long-range and static below Tgp (Fig. ||d). The zig-zag 
structure probably results from inter-layer couplings. 
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